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Abstract 

We investigate the role of exclusive channels in semi-inclusive electroproduction of pions and 
kaons. Using the QCD factorization theorem for hard exclusive processes we evaluate the 
cross sections for exclusive pseudoscalar and vector meson production in terms of generalized 
parton distributions and meson distribution amplitudes. We investigate the uncertainties 
arising from the modeling of the nonperturbative input quantities. Combining these results 
with available experimental data, we compare the cross sections for exclusive channels to 
that obtained from quark fragmentation in semi-inclusive deep inelastic scattering. We 
find that p° production is the only exclusive channel with significant contributions to semi- 
inclusive pion production at large z and moderate Q 2 . The corresponding contribution to 
kaon production from the decay of exclusively produced (ft and K* is rather small. 
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1 Introduction 



Electroproduction processes in the Bjorken regime probe the partonic structure of the nucleon. In- 
clusive deep inelastic scattering (DIS) provides extensive information about the sum of quark and 
antiquark distributions in the nucleon, and allows one to determine the gluon distribution from the 
observed scaling violations. More detailed information can be obtained from scattering experiments 
in which one or more hadrons are observed in the final state. There are two basic classes of such 
experiments. The first one is semi-inclusive deep inelastic scattering, in which one observes a single 
hadron, carrying a fraction z of the photon energy in the target rest frame, in a final state of large 
average multiplicity. A QCD factorization theorem states that the semi-inclusive cross section in the 
Bjorken limit can be expressed in terms of distribution functions for quarks, antiquarks and gluons 
in the target and of the corresponding fragmentation functions into the observed hadron. This allows 
one to tag the active parton via its fragmentation properties and has recently been used for a flavor 
decomposition of polarized quark and antiquark distributions in the semi-inclusive production of pi- 
ons and kaons [Q. In addition, measurements of azimuthal asymmetries in semi-inclusive pion and 
kaon production, such as the Collins and Sivers asymmetries for a transversely polarized target [2j, 
provide interesting information about the distribution of the spin and transverse momentum carried 
by quarks and antiquarks in the nucleon. 

The second class are exclusive scattering processes, in which the final state is a recoiling baryon 
B together with a single meson or a few-meson system carrying nearly the full photon energy in the 
target rest frame. A QCD factorization theorem states that in the Bjorken limit, and for longitudinal 
photon polarization, the amplitudes of such processes are calculable in terms of the light-cone distri- 
bution amplitude of the produced meson and generalized parton distributions (GPDs) for the p — > B 
transition. Generalized parton distributions provide a wealth of information on the parton structure 
of the nucleon, in particular about the spatial distribution of partons in the transverse plane and 
about quark orbital angular momentum, see e.g. the reviews [31 EJ 03 El Hj • 

The Bjorken limit for semi-inclusive electroproduction implies a large average multiplicity of the 
produced hadronic system. In practice, however, there can be situations in which individual exclusive 
channels play an important role. In fixed-target experiments the limited photon energy restricts 
the phase space for quark fragmentation, in particular at large z, and for relatively low photon 
virtuality Q 2 the suppression of individual exclusive channels due to the faster drop of the exclusive 
cross sections with Q 2 may not yet be effective. In particular, phenomenological studies suggest 
that a large contribution to tt^ 1 production comes from exclusive p° production, with subsequent 
decay p — > tt + tt~ [S1U5]. Fortunately, the cross sections for exclusive p° electroproduction has been 
measured by the HERMES and CLAS experiments, including the cross section ratio for longitudinal 
and transverse photons jlUl 1111 1121 1131 114j . It is natural to ask whether the strange vector mesons </> 
and K* play an equally important role in semi-inclusive kaon production and whether other exclusive 
channels may be important, too. A quantitative answer to these questions will help to delineate 
the limits of the kinematic region where semi-inclusive data can be analyzed using the factorization 
theorem. 

In this paper we investigate the role of exclusive channels in the semi-inclusive production of pions 
and kaons on the basis of the factorization theorem for hard exclusive processes. Our investigation 
consists of two parts. Firstly, we evaluate the longitudinal cross section for the exclusive production of 
pseudoscalar and vector mesons in the leading-twist approximation and at leading order in the strong 
coupling, focusing on tt, K and p, <fi, K*. We explore uncertainties in the obtained cross sections, 
in particular those due to the generalized parton distributions, which are still largely unknown and 
need to be modeled. Such uncertainties will persist if higher-order and higher-twist corrections are 
included. Seen from a different perspective, they indicate to which extent exclusive meson production 
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is sensitive to GPDs and thus interesting in its own right. It is known that exclusive meson production 
cross sections at moderate Q 2 are affected by substantial power corrections. For the production of 
p°, (p and 7r + there is data or preliminary data, which we will use to assess the quantitative validity 
of our calculated cross sections. Secondly, we evaluate the contribution of these exclusive channels to 
semi-inclusive production of ir and K, and compare them with the results obtained from leading-twist 
quark fragmentation. For the exclusive meson production cross sections we rely on experimental data 
where possible, and only use our leading-twist calculation to estimate the ratio of cross sections for 
measured and unmeasured channels. 

The paper is organized as follows. In Sects. GO and 01 we summarize the leading-twist description 
of exclusive meson production and describe the models for the GPDs used in our investigation. 
An analysis of pseudoscalar and vector meson production channels is then given in Sects. |1] and 
03 In Sect. El we briefly discuss the limitations of our leading-order results and compare them with 
experimental data. The contribution of exclusive channels to semi-inclusive meson production is 
compared with leading- twist quark fragmentation in Sect. and in Sect. 00 we summarize our main 
results. Some important technical details are given in two appendices. 



2 Exclusive meson production in the leading-twist approximation 

Let us consider the exclusive electroproduction process 

e(k) + p(p) - e(k') + M(q') + B(p'), (1) 

where M is a meson and B a baryon, and where four-momenta are indicated in parentheses. Through- 
out this work we assume beam and target to be unpolarized. We write q = k — k' , A = p' — p, and 
use the standard kinematic variables 



/ 



Q 1 



W z 



(P + Q? 



xb = Q 2 /{2pq), y = (pq)/(pk). 



(2) 



We respectively write m p , tub, tum for the masses of the proton, the baryon B, and the meson M. 
With Hand's convention JB] f° r the virtual photon flux, the electroproduction cross section is given 
by 



da(ep -> eMB) 



a. 



y z 1 - x B 1 



dQ 2 dxB dt 2tt 1 — 6 xb Q 2 

in terms of the cross sections do~T/dt and dai/dt of the -y*p - 
longitudinal 7*, where 

1 - V ~ {yxBm p /Q) 2 



da T da L ] 

> MB subprocess for transverse and 

(4) 



l-y + y 2 /2 + {yx B m p /Q) 2 

is the ratio of longitudinal to transverse photon flux. 

In the generalized Bjorken limit of large Q 2 at fixed xb and fixed t, the process amplitude factorizes 
into a hard-scattering kernel convoluted with generalized parton distributions for the p — > B transition 
and with the distribution amplitude of the meson ^S]. Example graphs are shown in Fig. ^ In this 
limit the longitudinal cross section can be written as 



da L 
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for vector mesons, and as 
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Figure 1: Example graphs for exclusive meson production in the generalized Bjorken limit. The large 
blob denotes the p — > B transition GPD and the small one the meson distribution amplitude. 

for pseudoscalar mesons M. The transverse cross section dax/dt is power suppressed by 1/Q 2 com- 
pared with da^jdt. Here we have in addition used the skewness variable 



x B 



(p-p')(q + q') 



(7) 



where the approximation holds in the generalized Bjorken limit. Note that the prefactor in © and 
© can be rewritten as x B /(l — x B ) = 4£ 2 /(l — £ 2 ). The quantities Tt, £ and H, £ are specific for 
each channel. Throughout this work we will take their leading order approximations in a s . To be 
specific, we have 
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(8) 



for 7*p — ► isT* + A and 7*p — > K + A, respectively, with analogous expressions for £k*+a an d £r+a- 
The GPDs for the p — ► A transition are defined as 
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m A + m p 



1 f dz~ 



ixP + z 



2tt 



(A| S -(-iz) 7 + 75 n(iz) 



z+=0,z T =0 



2P^ 



H p ^ A (x, £, t) u7 + 7 5 u + E p ^ A (x, £, t) u 



7 5 A+ 
mA + mr 



-u 



(9) 



where we use light-cone coordinates = (w° ± v 3 )/^/2 and ut = (f 1 ,^ 2 ) for a four-vector v and 
assume light-cone gauge A + = 0. For brevity we have not displayed the momentum and polarization 
dependence of the baryon spinors on the right-hand sides. GPDs for other transitions are defined in 
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Table 1: Combinations of proton GPDs to be used for various channels 7p — > MB at the place of 
2H p ^\(x, £, t) + H p ^\(—x, £, t) or 2H p ^\(x, £, t) + H p ^\(—x, £, t) in (jSJ). All distributions are to be 
evaluated at arguments x, £, t, with iif 9 , iif 9 and i/ 5 as defined in |S] and if', JH"* given above 1)121) . 



p°p ( [2F U + F d ] + [2£P + H d ] + | x" 1 ^ ) 

ujp 1 f [2iJ u - F d ] + [2iP - + | ar 1 ^ ) 

K*+A --^(2[2i/ u - F d - F s ] - [2iJ" - H d - if 5 ]) 

K*+Y,° -±(2[H d - H s ] - [H d - iP]) 

K*° £+ [H d - H s ] + [H 3 - H s ] 

c/)p -([H s + H s ] + \ x - l H9^j 

TT+n 2[H U - H d ] + [H a - H d ] 

TrOp -t= ( [2H U + H d ] - [2H n + H d ] ) 

A'+A -75 ( 2 i 2 H u - H d - H s ] + [2H n - H d - H~ s ] ) 

K+T, --^(2 [H d — H s ] + [H d — H s ]^j 

A°S + [H d - H s ] - [H d - H s ] 



full analogy. The integrals over meson distribution amplitudes in Q can be expressed as 



z[l-z) 



1 + °2n 



n=l 



through their coefficients in the expansion 
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z(l - z) 



6 ^2 a2n ~ 
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(10) 



z) = 62(1 - z) [l + jr a n Cl'\2z - 1) 



n=l 



(11) 



on Gegenbauer polynomials, where z is the light-cone momentum fraction carried by the quark in 
the meson. Note that odd Gegenbauer coefficients a,2n-i describe an asymmetry in the momentum 
distribution of the quark and antiquark in the meson. They can be nonzero for K and K* due 
to flavor SU(3) breaking. In (JHJ) to fj 1 1|1 we have not displayed the logarithmic dependence on the 
renormalization scale in a s and on the factorization scale in the GPDs and the distribution amplitudes. 

Using flavor SU(3) symmetry one can relate the transition GPDs from the proton to a hyperon 
to the distributions H q (x,^,t) for quark flavor q in the proton [3J Ej. This gives in particular 
H p ^ A = — [2H U -H d - H s ]/V~6 and an analo gous relation for H p ^\. We use these relations through- 
out, except for E, where there are large effects of SU(3) breaking as we shall see below. Results 
analogous to (jSJ) hold for all meson channels we consider, see e.g. jUElEl) an d we have collected 
the relevant combinations of GPDs in Tabled There we have introduced H q (x,^,t) = —H q (—x,^,t) 
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and H q (x,£,t) = H q {— t), so that for x > we have simple forward limits 

H q (x, 0, 0) = q(x), H q (x, 0, 0) = q(x), H q (x, 0, 0) = Aq(x), H q (x, 0, 0) = Aq(x) (12) 

in terms of the unpolarized and polarized quark and antiquark densities in the proton. For gluons 
we have H 9 (a;, 0,0) = xg(x), which is the origin of the additional factors x~ l in the entries for p°, 
lo, (j). In addition to the replacements in Table ^ one has of course to take the appropriate meson 
distribution amplitude and meson decay constants in JSJ. For the latter we will take f n = 131 MeV, 
f K = 160 MeV, and 

f p = 209 MeV, f u = 187 MeV, = 221 MeV, f K * = 218 MeV (13) 

from |18| . 

For a s in (jHJ) we will take the one-loop running coupling at the scale Q 2 , with three active quark 
flavors and Aqcd = 200 MeV. This gives a s = 0.34 at Q 2 = 2.5 GeV 2 , where we will show most of 
our numerical results. We will not attempt more refined choices of renormalization scale, as were for 
instance explored in ^U], since our principal use of the leading-order calculation will be to describe 
the relative size of cross sections for different exclusive channels. 



3 Modeling the generalized parton distributions 

For the calculation of exclusive cross sections we use simple models of GPDs. They have been 
developed in |20| 0] and been used in most phenomenological analyses so far. Our aim here is not 
to improve on these models, but instead to see by how much predictions can vary within the given 
framework. We take a factorizing t dependence for H and H, 

H q (x,t,t) = H q (x,OG A (t)/G A (0), (14) 

where Ff(t) is the electromagnetic Dirac form factor of the proton and GA{t) the isovector axial form 
factor of the nucleon. A more refined version of the model would take different combinations of the 
proton and neutron form factors for H u and H d , but for the low values of t dominating integrated 
cross sections, F™(t) is much smaller than Ff(t) and we simply neglect it. In this sense (|14|) is 
consistent with the sum rule for the first moment / dx H q (x,^,t). The ansatz for H q is consistent 
with the sum rule for J dx H q (x,^,t) to the extent that the (unknown) isoscalar axial form factor 
has the same t dependence as the isovector one. In our numerical evaluations we take the familiar 
parameterizations 



Am 2 p - 2M 1 GAt)_ x 

lU 4m2-t [l-t/(0.71GeV 2 )] 2 ' G A (0) [1 - 1/(1.05 GeV 2 )] 2 



We note that for the gluon distribution H 9 there is no reason a priori to take the electromagnetic form 
factor Ff(t) in the ansatz JHJ. It turns out, however, that Ff(t) is well approximated by a dipole 
form Ff (t) = [1 - i/(0.98 GeV 2 )] -2 for t up to about 3 GeV 2 [2T] and thus close to the two-gluon form 
factor advocated in [22] • 

It is rather certain that the ansatz (|14[) is too simple and can at best reflect the correct t dependence 
in a limited range of x and £ [U 1211 I23j . For x and £ in the valence region, say above 0.2, the 
decrease of GPDs with t is most likely less steep than the one of Ff (t) and G A (t). Whereas there arc 
phenomenological constraints of the t behavior of valence quark GPDs |2J and for gluons at small x 
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|22| . the behavior for sea quarks is largely unknown, and sea quarks are important for the xb region 
around 0.1 we will be mostly concerned with. Furthermore, the t dependence of meson production at 
moderate Q 2 is strongly affected by power corrections, as is for instance seen in the Q 2 dependence of 
the logarithmic slope B = (d/dt) log(da / dt)\t=o for p° production at very high energies We note 
that our ansatz (j!4| gives a slope parameter B ~ 4GeV 2 , which may be quite realistic for xb around 
0.1. Furthermore, cross section ratios should be less affected by the insufficiency of our ansatz, since 
they are sensitive only to the relative t dependence of the contributions from different quark flavors 
and from gluons. 

For the t independent functions in (|T4"|) we use the double distribution based ansatz of [20], whose 
ingredients are the usual parton densities at a given factorization scale fi and a so-called profile 
parameter 6, where /i and b are to be regarded as free parameters of the model. Explicit expressions 
are given in App. EI We will not take into account the evolution of GPDs, which should not be too 
problematic since our numerical applications will stay within a rather narrow range of Q 2 . 

The modeling of the nucleon helicity-flip distributions E q and E 9 is still at an early stage of 
development, with the most advanced considerations focused on the valence quark domain [H I21|. 
Fortunately, contributions from E enter the unpolarized meson production cross section © with pref- 
actors that are quite small in the kinematics we are most interested in. Following the argumentation 
of |25j that E is not much larger than H for a given parton species, we hence neglect E altogether in 
our cross section estimates. 

The distributions E cannot be neglected since they receive contributions proportional to that 
compensate the kinematic prefactors in the cross section (jHJ). We model them as 



E p ^ A (x,£,t) 



0(\x\<£) fx + f\ 2m p f n g nNN A 



2£ 
2£ 



2£ 

x + j 

2£ 



m 2 n -t A 2 -t ' 
(m p + m A )f K g KNA A 2 



m 



K 



in 



K 



t 



A 2 -t 



(16) 



where the distribution amplitudes (ft are the same as those introduced above. For the coupling 
constants we take the value g n NN — 277i p Ga(0)/ f n ~ 14.7 from the Goldberger-Treiman relation and 
9k n a ~ —13.3 from Continued to the points t = m 2 or t = m\ in the unphysical region, the 
expressions (|16|) become the well-known results from pion or kaon exchange |261 1271 05] . These can 
only be expected to be good approximations for t close to the squared meson masses, and for —t of 
several 0.1 GeV 2 are to be regarded as extrapolations. In (|16|) we have included form factors that cut 
off the l/(m 2 M — t) behavior of the pure pole terms when —t becomes large. As default value for the 
cut-off mass we will take A = 0.8 GeV and study the sensitivity of results to the precise value of 
this parameter. We note that for A = 0.6 GeV and — t < 1 GeV 2 the above form of E p _> n differs by 
less than 15% from the corresponding term calculated in the chiral quark-soliton model as given in 
Eq. (4.39) of EHI 

With this model for E the longitudinal cross section for ^y*p — ► ir + n takes the form 



d(T L 
dt 



x B 



Q b 1 — xb 



(1 - e)\m, t)\ 2 - 2m v i Re7Y(£, t) Q 2 F n (Q 2 
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(17) 
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Figure 2: Leading-twist cross sections for ^)* L p — > 7r + n (left) and ^* L p — ► K + A (right) for Q 2 = 
2.5 GeV 2 . An upper cut on —t of 1 GeV 2 has been applied here and in all further plots of this paper. 
Shown are the individual contributions from \7i\ 2 and \£\ 2 , and their coherent sum according to (JHJ). 



is the electromagnetic pion form factor to leading order in a s and 1/Q 2 . An similar expression 
involving F K +(Q 2 ) is obtained for j*p — > K + A according to © and Q16JI . Note that the the \£\ 2 term 
in duL/dt has no xb dependence other than from the explicit factor x B /(l — xb)- Within our model 
the \H\ 2 term reflects the behavior of the polarized parton distributions at momentum fractions of 
order £, and its contribution to da^jdi can very roughly be represented by [^Ag(^)] 2 . 



4 Exclusive pseudoscalar meson production 

In this and the next section we present numerical results for cross sections of exclusive meson pro- 
duction. Our main focus is to compare the rates for different production channels and to investigate 
model uncertainties. Comparison with data in Sect. will allow us to estimate the shortcomings of 
the leading approximation in 1/Q 2 and in a s , on which our calculations are based. 

The factorization theorem for exclusive meson production requires t to be much smaller than Q 2 . 
For definiteness we will give all meson cross sections in this paper integrated over —t from its smallest 
kinematically allowed value —to to an upper limit of lGeV 2 . In generalized Bjorken kinematics we 
have 

2g 2 (m|+m 2 ) + 2^(m|-m 2 ) 
- t « ^— ^ (19) 

with £ defined in (|7|). For low enough xb most of the cross section should be accumulated in this t 
region, whereas for large xb our cross sections decrease to the extent that —to approaches lGeV 2 . 

To begin with, let us investigate the relative importance of the contributions from Ti and £ to 
7r + and K + production with our model assumptions. As is seen in Fig. [2J exclusive tt + production 
receives comparable contributions from both the \Ti.\ 2 and the \£\ 2 term in Q, as well as from the 
interference term proportional to Jie(£*7i). Note that the relative weight of the contributions is 
different for dai /dt, where it strongly depends on t given the characteristic t dependence of the pion 
pole term (j!6|) . In K + production the influence of £ is less prominent, since the pole factor (m 2 K — t)~ l 
gives much less enhancement at small t than (to 2 — t) . 

To obtain the curves in Fig. [21 we have made a number of choices in the nonperturbative input 
to the cross section, which we now discuss in turn. In Fig. |31we show how the cross section changes 
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Figure 3: Leading-twist cross section for j^p — > ir + n (left) and "i\p — > K + K (right) at Q 2 = 2.5 GeV 2 
obtained with different values of the parameter A in the form factor multiplying the pion or kaon pole 
contribution in 



when we vary the parameter A in our model for E, where A = 1.3 GeV represents an upper limit 
of the values discussed in the phenomenological study j^Hj, and A = 0.6 GeV approximates the form 
factor dependence obtained for the pion pole contribution in [2H1, as discussed at the end of Sect. |21 
Omitting the form factor altogether (tantamount to setting A — ► oo) the 7r + cross section would 
increase by more than a factor 1.4 and the K + cross section by more than a factor 1.7 compared with 
our default choice A = 0.8 GeV. Also, the cross sections would considerably increase when raising the 
upper cutoff in the —t integration above 1 GeV 2 . In other words, the cross section would then receive 
substantial contributions from values of t far away from the region where the pion or kaon pole term 
can be regarded reasonable approximation of E. 

For the pion distribution amplitude we have taken the asymptotic form <f> n (z) = 6z(l — z) under 
scale evolution, which is close to what can be extracted from data on the 7~7r transition form factor, 
see e.g. [313 EI] The study in [3^ quotes limits on 02 + 04 at scale \x = 1 GeV of about ±0.3 if all 
other Gegenbauer coefficients are set to zero. This corresponds to a change of the 7*p — > ir + n cross 
section by a factor (1 + 02 + a^) 2 between 0.5 and 1.7. For the K + distribution amplitude we have 
taken the asymptotic form as well. Figure 0] shows how the K + cross section changes if instead one 
takes a\ = —0.05 and 02 = 0.27 at \x = 1 GeV from the QCD sum rule calculation [221 • This value of 
a\ is compatible with the findings of |33| . 

For our model of H we have taken a double distribution ansatz with a profile parameter 6 = 1 
(see Sect. |21 and App.^J. Taking 6 = 2 instead would decrease the K + cross section by a factor of 
approximately 0.6. The pion cross section changes less, because of the relative weight of Ti and E. A 
more important source of uncertainty is however due to the polarized quark densities used as input 
to model H. As a default we have used the LO parameterization from (21] at a scale ji = 1 GeV. 
Using instead the LO parameterization in scenario 1 of [22] at the same scale, the K + cross section 
changes as shown in Fig. 21 Note that any uncertainty on parton distributions is amplified in the 
meson production cross section, where GPDs appear squared. 

Let us now comment on other pseudoscalar channels. The cross sections for ^*p — > K + TP is 
about an order of magnitude smaller than the one for 7*2? — ► K + A, as is seen in the numerical study 
of |26| . For the contribution from 7i this can be understood from the flavor structure in Tabled 
where for a rough estimate one may concentrate on the dominant terms H u and H d . For current 
parameterizations of polarized parton densities the combination [2H U — H d ]/\^3 for A production is 
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Figure 4: Left: Leading-twist cross section for ^p — > K + K at Q 2 = 2.5 GeV 2 calculated with the 
asymptotic kaon distribution amplitude and with the one from Khodjamirian, Mannel and Melcher 
(KMM) [321 • Right: The same cross section calculated with different parton distributions in the 
model for H. The distributions from Leader, Sidorov and Stamenov (LSS) |34j and from Bliimlein 
and Bottcher (BB) |3S] are taken at a scale (i = 1 GeV. 

clearly larger than H d in the analogous expression for the E° channel. Concerning the contribution 
from £ , the coupling gxNT, is about three times smaller than gxNA appearing in (|16|). see HHI- Along 
the same lines one can see that the cross section for j*p — > K°Ti + is of similar size as the one for 

The channel 7*p — ► ir°p does not receive contributions from the pion pole term in (|16[) because of 
charge conjugation invariance, so that in our model it is entirely given by the contribution from 7i. In 
Table^we see that the combination 2H U + H d for 7r° production is to be compared with \^2[H U — H d ] 
for j*p — ► 7r + n, which is of comparable size. One thus expects the 7r° cross section to be similar to 
the \7i\ 2 part of the 7r + cross section. 

The exclusive channels j*p — > rjp and j*p — > rjp involve the combination 2H U — H d instead 
of 2H U + H d in the 7r° case, which is somewhat larger because the polarized distributions Ait and 
Ad have opposite sign. The strange quark contribution to these channels involves H s — H s , which 
vanishes in our model with polarized parton distributions satisfying As(x) = As(x). A quantitative 
analysis requires the appropriate decay constants for the r\ and rf , see for instance [301; but one can 
in general expect comparable cross sections for the 7r°, r\ and rj' channels. 

5 Exclusive vector meson production 

Within our model for the GPDs, the cross section for vector meson production is sensitive to unpo- 
larized quark and gluon densities. To obtain an indication of uncertainties we have compared results 
with the LO distributions from CTEQ6 £21 and the LO distributions from MRST2001 [SHI - For 
consistency we need LO rather than NLO parton densities, which unfortunately are not available in 
several of the most recent parton fits. We have checked that the NLO distributions from MRST2001 
are in good agreement with those in the MRST2002 and MRST2004 analyses |39j for quark and 
antiquark densities down to about x ~ 1CT 2 and for the gluon density down to about x ~ 1CP 1 . 
Comparing the LO distributions of CTEQ6 and MRST2001 at a scale /i 2 = 1.2 GeV 2 (which is the 
lowest value accepted by the code for the MRST2001 parameterization) we find that the CTEQ6 
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Figure 5: Leading-twist cross section for — > p°p at Q 2 = 2.5 GeV 2 . Shown are the individual 
contributions from quark and gluon distributions and their coherent sum. The parton densities in 
the double distribution model are taken at scale fi 2 = 1.2 GeV 2 for the upper and at p 2 = 2.5 GeV 2 
for the lower plots. 



gluon is larger for x< 10 and smaller for x> 10 . The u quark distribution is quite similar in the 
two sets for x > 10 -2 , whereas the s quark is significantly smaller for CTEQ6. The distributions for 

d, u, d are quite similar for x> 10 -1 and larger for CTEQ6 at smaller x. The LO parameterization 
of Alekhin [3U] has significantly larger u and u distributions and a smaller gluon than the two other 
sets. At p? = 1.2 GeV 2 it has however almost no strange quarks in the proton, which we do not 
consider physically very plausible and which is in clear contrast with the results of CTEQ6, where a 
dedicated analysis of data constraining the strangeness distribution was performed. Since our study 
is crucially dependent on the flavor structure of parton distributions, we have therefore not used |4Uj . 
Comparing the different parton sets at the higher scale p? = 2.5 GeV 2 we find a very similar picture. 

In the double distribution model of GPDs we take the profile parameter 6 = 2 for both quark 
and gluon distributions. For all mesons we take the asymptotic shape of the distribution amplitude, 
given that no direct experimental information is available for them. Theoretical estimates do not give 
stronger deviations from the asymptotic form than those we discussed for pseudoscalar mesons, see 

e. g. the compilation in [ISj. In Fig. Owe show the individual contributions from quark and gluon 
distributions to the p° cross section as well as their coherent sum. The clear difference between the 
CTEQ6 and MRST2001 result reflects the current uncertainty on the usual parton densities at low 



12 



0.1 0.2 0.3 0.4 0.5 0.6 

X B 



0.1 0.2 0.3 0.4 0.5 0.6 



Figure 6: Leading-twist cross sections for ^y L p — > p°p and "y^p — > up (left) and for ^y L p — > p + n 
(right) at Q 2 = 2.5 GeV 2 , obtained with the MRST2001 parton densities taken at p 2 = 1.2 GeV 2 . 

scales in the relevant range of x. A striking feature is the clear dominance of the gluon distribution up 
to quite high values of xb- Note that with our model of GPDs the convolutions Ti. q (^,t) and 7i 9 {^,t) 
are sensitive to the forward parton distributions in a certain range of momentum fractions around £ 
(see App.^). The strong dominance of gluon over quark distributions at small momentum fractions 
thus still shows its effect at xb values above 10 -1 . Note that we have taken the same t dependence 
for quark and gluon GPDs in our model (|14j) . lacking phenomenological evidence to the contrary. 
Comparison of the t dependence e.g. for p° and p + production in equal kinematics could be of help 
here. 

In our numerical calculations we have calculated the integrals ()36|) and (|37|) for the meson produc- 
tion amplitude with a lower cutoff at momentum fractions x = 10 -4 . The cross section for xb = 0.05 
changes by less than 2% if we cut off at 10 -5 or at 10~ 3 . It is diminished by about 10% with a 
cutoff at 0.005, which gives an indication of the relevance of momentum fractions which are an order 
of magnitude smaller than the xb of the process in this model. Similar changes are observed for 
the individual quark and gluon distributions. We note that if we take a profile parameter 6=1 for 
quarks, the quark contribution to the cross section at xb = 0.05 decreases by 10% when moving the 
cutoff on x from x = 10~ 4 to x = 10~ 3 and by 35% when moving it from x = 10~ 4 to 0.005. Such 
a strong dependence on momentum fractions well below xb seems difficult to understand in physical 
terms. We note that in the sea quark sector there are no strong theoretical arguments for taking 
6=1, see Sect. 4.4 of 0. 

In Fig. [21 we also observe a significant change of the gluon contribution to the cross section when 
changing the scale of the parton distributions in the double distribution ansatz (|32|) . In contrast, the 
quark contribution changes by at most a factor of 1.3, reflecting the relatively weak scale evolution of 
quark and antiquark distributions compared with gluons in the relevant kinematic region. Changing 
the scale of the forward distributions in the double distribution model (j3*2*|) gives a rough indication 
of how the actual GPDs evolve with p 2 [20], so that the strong increase with p 2 seen in Fig. \5\ 
reflects a strong scale uncertainty of the leading-order approximation in a s for channels involving 
gluon exchange. A full NLO analysis of meson production is possible using the results of [U but 
beyond the scope of this work. We will use the smaller scale p 2 = 1.2 GeV 2 in our further studies, 
because the internal virtualities in the hard-scattering graphs of Fig. n are typically smaller than Q 2 
(see also the study |46j of relevant scales in the small- x limit). Furthermore, the MRST2001 set gives 
a better description for the ratio of (fr and p° production with our model (see below) and will hence 
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be our default choice in the following. 

In Fig. Elwe show the production cross sections for p + , p° and oj. The xb behavior of the p + 
cross section roughly follows the one of £ 2 [it(£) — d(£) + u(£) — d(^)] 2 , which is a flavor nonsinglet 
combination and hence does not display the strong rise of sea quarks or gluons at small x. The clear 
suppression of u production compared with the p° is a consequence of the relative factor in the gluon 
contribution (see Table ^) and at large xb of the relative size of the flavor combination 2H U — H d 
compared with 2H U + H d . We remark that the exclusive channel 7*j> — > j^V a l so contributes to 
semi-inclusive production of ir + , tt" and vr°. It involves the combination 2H U + H d - [2H U + H d ], 
where sea quarks drop out, so that one may expect a cross section of similar size as for p + production. 
A numerical estimate would however require knowledge of quark and gluon distribution amplitudes 
of the /2, see |42| . and is beyond the scope of this work. 

Figure [7] shows our results for K* + and K*° production. In contrast to K + production, the cross 
section for the A channel is not much larger than for the S° channel. Consulting Tabled we see that 
this is because the contributions from u and d quarks partially cancel in 2H U — H d whereas they add 
in 2H U — H d . We remark that the results for K* + and K*° production decrease by less than 25% 
when instead of MRST2001 we take the CTEQ6 parameterization. The uncertainty due to knowledge 
of the parton distributions is hence much less than for the gluon dominated channels. 

Results for (ft production are shown in Fig. |H1 The dominance of the gluon over the strange quark 
contribution is clearly seen, although strange quarks are not entirely negligible with the MRST2001 
parameterization. 1 Since gluons dominate for most xb, we see the same trend concerning differences 
between the parameterizations and the choice of scale as for p° production. The ratio of <jl for (ft 
and p° production is shown in Fig. El where the dependence on p 2 is seen to be much milder since it 
partially cancels in the ratio. The difference between CTEQ6 and MRST2001 is still significant and 
mainly due to the difference in the gluon distributions. 

Preliminary data from HERMES [111 ILSj give a ratio of about 0.08 for the cross sections of (ft and 
p° production for x B = 0.09 and Q 2 = 2.46 GeV 2 and for x B = 0.13 and Q 2 = 3.5 GeV 2 . These data 
contain a significant contribution from cry, and preliminary HERMES data |11 | I12 | IT3"] suggest that 
R = ctl/vt may be slightly smaller for (ft than for p° production at the same Q 2 . The (ft to p° ratio 
for gl would then be somewhat larger than 0.08. In addition, one can expect that a narrower shape 

1 In the study 1251 strange quarks were neglected based on inspection of the CTEQ6 parameterization. 
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Figure 8: Leading-twist cross section for 7£p — > 4>p at Q 2 = 2.5 GeV 2 . Shown are the individual 
contributions from gluons and the coherent sum of gluons and strange quarks. The upper plots are 
for parton densities taken at scale /z 2 = 1.2 GeV 2 in the double distribution model, and the lower 
plots for parton densities taken at /i 2 = 2.5 GeV 2 . 



of the distribution amplitude and power corrections due to the strange quark mass would decrease 
the estimates in Fig. |§] 25 . 

A complete representation of GPDs includes in addition to the double distribution the so-called 
D-term . It vanishes in the forward limit £ = and does not affect the double distribution ansatz 
we are using. Its contribution to the GPDs can be expanded in Gegenbauer polynomials as 

(2 \ oo 
l-^)E4n +1 (*)^ 2 +1 (|), 

HUx,a,t) = 0(w<e)f (i-|H fx+i(*)<^ 3 (f) (20) 

for £ > 0. Such terms contribute to the real part of the convolution integrals needed in the meson 
production amplitudes as 

_ /■ d Hi(*, ( ,t) _ r dx f M _ 2 g 

D J-i i-x-ie J-i i-x-ie 
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Figure 9: Ratio of leading- twist longitudinal cross sections for <f> and p° production, obtained with 
different parton distributions taken at /i 2 = 1.2 GeV 2 (left) and with the MRST2001 distribution 
taken at different fi 2 (right). 

' n=0 

These terms give a £ independent contribution to TC(t;, t), in contrast to the contributions from the 
double distribution part, which very roughly follow the behavior of £(/(£), £<?(£) or £<?(£) and hence 
grow as £ becomes smaller. In |44j the first three coefficients in the quark D-term at t = have been 
extracted from a calculation in the chiral quark-soliton model of the nucleon, giving df(0) ~ —2.0, 
cfij(O) w —0.6, dg(0) « —0.2 and equal values for d quarks, referring to a scale fj, of a few GeV [4*5] . 
The gluon D-term is parametrically subleading at the low scale intrinsic to the model, but evolution 
to [i of a few GeV can give values similar in size to those for quarks. The values Ijj = if) = —5.6 
turn out to be similar in size and opposite in sign to the real parts of the corresponding integrals 
from the double distribution part in our model. The effect of such a D-term is however much weaker 
on the square \H\ 2 appearing the cross section, which is dominated by the imaginary parts of the 
integrals in a large region of xb- Taking the above values for the quark D-teim and I 9 D = —11.2 as 
an order-of-magnitude estimate, we find that the change of the various vector meson cross sections 
is at the 10% level for xq = 0.1 and not more than a factor 1.5 in either direction for xb below 0.3, 
where for definiteness we have taken the MRST2001 distributions at fi 2 = 1.2 GeV 2 . 

6 Comparison with data and discussion of power corrections 

Our calculation of meson production cross sections is based on the leading-twist approximation. It is 
known that corrections in l/Q 2 to leading-twist meson cross sections can be substantial for Q 2 of a 
few GeV 2 . A systematic treatment of such power corrections remains an unsolved problem. There 
is however a number of approaches that allow one to model particular sources of power corrections, 
see e.g. [5] for a discussion and references. For vector meson production, a considerable suppression 
of the leading-twist result at moderate Q 2 is found when including in the hard-scattering kernel 
the transverse momentum of the quarks in the meson |461 1471 148j . This means that the transverse 
resolution power of the virtual photon cannot be neglected compared with the transverse size of the 
meson. Similarly, the finite transverse momentum of the partons coming from the proton gives rise 
to power corrections, when it is included in the hard-scattering kernel. Estimating both effects in 
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Table 2: Experimental values cr^ xp of the longitudinal cross section for the production of p° ^U] and 
4> (111 1131 149j from HERMES, and the ratio between our leading- twist results a^ y (obtained with 
the MRST2001 distributions taken at p 2 = 1.2 GeV 2 ) and the data. The data for <j) production is 
preliminary. 



Q 2 [GeV 2 ] 


7*p — > p°p 
xb <r™ p [/ib] 




1*P- 
Q 2 [GeV 2 ] 


~* (preliminary data) 

exp r l i thy / exp 

x B cr L y [TLb\ a L y /a L p 


2.3 


0.1 0.21 ±0.04 


7.1 


2.3 


0.087 15.6 ±3.1 14.9 


2.3 


0.075 0.21 ± 0.04 


7.6 


3.8 


0.136 6.2 ± 1.24 5.5 


4.0 


0.16 0.09 ±0.02 


2.1 






4.0 


0.12 0.06 ±0.01 


3.5 







a calculation considering only quark GPDs [I7|; a suppression of the leading- twist cross section for 
p° production by factors of 3.3, 4.9 and 9.2 was found at Q 2 = 5 GeV 2 for xb = 0.3, 0.45 and 0.6, 
respectively. The recent study [IS] for small xb, where only gluon GPDs were retained and only the 
transverse quark momentum in the meson was taken into account, found corresponding suppression 
factors of 4.6 and 6.6 for respective values of xb = 2.95 x 10~ 3 and 10~ 4 at Q 2 = 4.8 GeV 2 . For 
Q 2 = 10.9 GeV 2 and xb = 4.3 x 10~ 3 this factor decreases to 1.9. The discrepancy of the calculation 
including power suppression and experimental data is less than 35% in all three cases. 2 

In Table |2 we compare our leading-oder results for p° and (j) production with data from HERMES 
|101 lll| I13| I49j . The discrepancy between our calculation and the p° data is well in the range of 
what can be explained by suppression from quark transverse momentum (considering in addition the 
uncertainties of our model for the GPDs). The stronger discrepancy with the 4> data corresponds to 
our overestimating the 4> to p° production ratio, discussed in the previous subsection. 

The CLAS collaboration has published results for p° production at xb = 0.31 and at xb = 0.38, 
with Q 2 values between 1.5 and 2.3 GeV 2 |14| . and for <fi production at xb = 0.29 and Q 2 = 1.7 GeV 2 
|50j . We consider that this kinematics, where the hadronic invariant mass W is below 2.3 GeV, is too 
close to threshold for comparison with a leading-twist calculation. A recent CLAS measurement 
of uj production at higher energies, with Q 2 up to 5.1 GeV 2 and W up to 2.8 GeV, found that helicity 
conservation between the 7* and the uj is strongly violated, in contrast with the predicted behavior 
in the large-Q 2 limit. This prevented the extraction of &l and was ascribed to a strong contribution 
from 7T° exchange (which is absent in the p° and (j) channels). 

Let us now turn to ir + production, where the situation is quite different. For the contribution 
from 7i one expects a similar suppression from quark transverse momentum as in the case of 7i in 
vector meson production, which was indeed found in the numerical study |47j . The pion exchange 
contribution from £ is described in terms of the pion form factor according to (|17|) , and this relation 
persists beyond the leading approximation in 1/Q 2 to the extent that the pion emitted from the 
nucleon is not too far from off-shell. The power corrections for £ are then the same as those for the 
pion form factor. The leading- twist expression (HEP for F^Q 2 ), with our choice of a s specified after 
(U2J), undershoots the data of [S3 by a factor 0.53 at Q 2 = 1 GeV 2 and a factor 0.41 at Q 2 = 1.6 GeV 2 . 

2 The leading-order formula (90) in |48) with which we obtained the suppression factors just quoted contains in addition 
an approximation for small xb, which should however not be the dominant effect comparing to the full calculation. 
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For Q 2 = 3.3 GeV 2 we find this factor to be between 0.34 and 0.77 within the large error bars in |53| . 
We will not attempt here to summarize the detailed theoretical and phenomenological work on the 
pion form factor, but remark that in addition to the leading-twist perturbative contribution there is 
a contribution from the Feynman mechanism, where the photon hits a quark carrying almost all of 
the pion momentum. This leads to a considerable enhancement over the leading-twist approximation. 
The calculations of F^iQ 2 ) in [21], which take this effect into account using different methods, give for 
instance results larger than our leading- twist value by factors between 2 and 4, even at Q 2 = 10 GeV 2 . 
Note that these factors are to be squared in the contribution of \£\ 2 to the ir + production cross section. 
For the production of the neutral pseudoscalars 7r°, 77, rf , where there is no pion exchange contribution, 
one expects that power corrections will decrease the cross section, similarly to the case of vector meson 
production. 

We have further compared our leading- twist calculation of eox for 7*p — > ir + n with preliminary 
data on ot + col from HERMES [55,. The HERMES data are presented for three different bins in 
xb, with the average values of Q 2 and xb for individual data points ranging from 1.5 GeV 2 and 0.1 
to 4.2 GeV 2 and 0.17 in the first bin, from 2.5 GeV 2 and 0.21 to 6.3 GeV 2 and 0.25 in the second bin, 
and from 4.5 GeV 2 and 0.34 to 10.5 GeV 2 and 0.45 in the third bin |56j . Averaging the ratio between 
theoretical and experimental cross sections for the data points in each bin, we find ratios of 0.42, 
0.19 and 0.12 in the first, second and third bin, respectively. The large discrepancy at large Q 2 and 
large xb (which are strongly correlated in the data) is striking, but not too surprising given the size 
of corrections just estimated for the pion form factor. The much better agreement at smaller Q 2 and 
xb might be accidental, given that we expect comparable contributions from H and £, for with the 
power corrections go in different directions. Help in clarifying this issue could come from the spin 
asymmetry for transverse target polarization, which gives access to the relative size of Ti and £ 

The case of tt + production (and also the findings in the u> channel mentioned above) show that 
there are specific power corrections which will not cancel in cross section ratios for different processes. 
The situation is however better for channels that are sufficiently similar, as the example of p° and 
4> production shows. Corrections due to quark transverse momentum (as well as the overall normal- 
ization uncertainty from the scale of a s in our leading-order calculation) will tend to cancel in that 
case. We hence expect that the overall pattern of differences between various meson cross sections 
we estimated at leading order will not be overturned in a more realistic treatment, given that these 
differences are largely controlled by the relevant combinations of quark and gluon distributions. 



7 Exclusive channels in semi-inclusive pion and kaon production 

In semi-inclusive hadron production one considers processes of the type 

e{k) + p(p) - e(k') + h(q h ) + X, (22) 

where h is a specified hadron and X an unspecified inclusive final state. A basic observable is the 
distribution of the produced hadron over the variable 

QhP , ON 

z = , (23) 

qp 

which measures the fraction of the virtual photon energy carried by the produced hadron in the target 
rest frame. In the Bjorken limit of large Q 2 at fixed xb and fixed z, semi-inclusive hadron production 
can be treated within a QCD factorization approach. The differential cross section factorizes into 
the distribution of partons of type i in the target, the cross section for the virtual photon scattering 
off this parton, and the fragmentation function D % ^ h {z) describing the fragmentation of the struck 
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Figure 10: Contributions of exclusive channels to semi-inclusive pion and kaon production, calculated 
at leading order, (a) Direct exclusive production of pseudoscalar mesons, (b) Exclusive production 
of a vector meson with subsequent decay into pseudoscalar s. 

parton into the hadron h, which carries a fraction z of its longitudinal momentum. To leading order 
in a s one has 

da n y 2 1 



2vra em J- — Y, 4 \q{x B )D^\z) + q{x B )D^\z)] , (24) 



dQ 2 dxB dz 1 — e xbQ 4 

where the sum is over quark flavors. Note that 1)24(1 has the same Q 2 -dependence as the inclusive 
DIS cross section in the Bjorken regime. The fragmentation functions are process-independent and 
describe not only semi-inclusive DIS but also e + e~ annihilation into hadrons and the distribution of 
leading hadrons in high-p^ jets. Their scale evolution is governed by evolution equations analogous 
to the DGLAP equations for the parton distribution functions. Various parameterizations of the 
fragmentation functions have been presented in the literature, which fit e + e~ annihilation and semi- 
inclusive DIS data at higher scales. 

The derivation of semi-inclusive factorization relies on the fact that in the Bjorken limit the 
inclusive final state X has a large invariant mass 

m 2 x = Q 2 ^—^+ml-(q-q h )\ (25) 

XB 

and thus a large average multiplicity (note that the squared momentum transfer (q — qh) 2 is always 
negative). The semi-inclusive cross section is thus obtained by summing over many individual chan- 
nels. In practice m 2 x is not very large for Q 2 values of a few GeV 2 typical of fixed-target experiments, 
for instance at Jefferson Lab or HERMES, especially at high z. At the same time, for moderate Q 2 
the suppression of the cross sections for exclusive channels relative to the semi-inclusive cross section 
(see below) may not yet be effective. One thus can reach a situation in which the cross sections 
of individual exclusive channels becomes comparable to the semi-inclusive one. It is interesting to 
compare the semi-inclusive cross section (|24j) with the cross sections of exclusive channels contribut- 
ing to semi-inclusive production. In the following we investigate the role of exclusive channels in 
semi-inclusive ir and K production on a proton target. We study two types of exclusive channels (see 
Fig. 110 j) at a quantitative level: 

i) direct exclusive production of pseudoscalar mesons, j*p — > ir + n and 7*7? — ► K + A, 

ii ) exclusive production of neutral or charged vector mesons p, (ft, K* with subsequent decay into 
pseudoscalars. 

In the direct production of pseudoscalar mesons (Fig. llUb ). the energy fraction z carried by the 
produced meson is related to the invariant momentum transfer to the nucleon t by 

mo — m 2 — t 

l-z = x B 3 Q2 P (26) 
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according to (|25jl. Exclusive production in the limit of large Q 2 at fixed xb and fixed t thus corresponds 
to values of z very close to 1. For example, tt + and K + production corresponds to z > 0.94 in 
typical HERMES kinematics of x B = 0.1 and Q 2 = 2.5 GeV 2 with a maximum momentum transfer 
\t\ = 1 GeV 2 (since the cross section drops rapidly with t, most events have z values yet closer to 
unity). Such exclusive contributions can usually be separated experimentally from the semi- inclusive 
events at smaller values of z. 

The situation is different for the contribution to semi-inclusive production resulting from the decay 
of exclusively produced vector mesons (Fig. HOb ). Since the decay products share the energy of the 
vector meson, such contributions result in an extended z distribution for the pion or kaon, even in the 
Bjorken limit. With the approximations described in Appendix[Bl the z spectrum of the pseudoscalar 
meson Pi from the decay V — > P1P2 can be written as 

da(ep -> Pi + P 2 B) a em y 2 1 - xb 



dQ 2 dxB dz 2ir 1 — e xbQ 2 

with 



ecTL^p VB) D L (z) + a T (j*p VB) D T {z) (27) 



= ^3 (* ~ z o) 2 , D T {z) = A (z - z 1 )(z 2 - z). (28) 



Here z\ < z < z 2 with 



zo = — — , z\ = z - C , z 2 = z + ( , C = (29) 

my my 

up to corrections of order XBm 2 /Q 2 . For brevity we have not explicitly indicated the dependence of 
Dl and on xb and Q 2 due to these corrections. The energy and three-momentum of Pi in the 
rest frame of the vector meson 

E = my + mpi ~ m 2 P2 
2my 



im v + mp l + m P2 — 2 {m v m P1 + m v mp 2 + m P1 mp 2 ) 
^ = ~ 2^ (30) 

depend only on the meson masses, and so do zq, z\ and z 2 in the limit of large Q 2 . In particular, the 
smallest and largest possible values of z for pions from the decay p — ► tttt are z\ = 0.04 and z 2 = 0.96 
in that limit. The corresponding values for kaons from <j) — > KK are Z\ = 0.37 and z 2 = 0.63. For 
the kaon from K* decay one has z\ = 0.32 and z 2 = 0.96, and for the pion from K* decay one has 
zi = 0.04 and z 2 = 0.68. 

According to (|5 J ) and (|27|) . the contribution of exclusive vector meson production to the cross 
section da / (dQ 2 dxs dz) asymptotically scales as l/Q s at fixed xb and z and is thus suppressed by 
l/Q 4 compared with the leading behavior (|24j) of the semi-inclusive cross section. Notice that (|24|) 
corresponds to transverse photon polarization, with contributions to the longitudinal cross arising at 
the level of a s and of 1/Q 2 corrections, just as in the familiar case of inclusive DIS. The situation 
is opposite for hard exclusive meson production, where <tl dominates over ax in the large-Q 2 limit. 
Measurements show however that at Q 2 in the few GeV 2 region the ratio R = olI&t is still of order 
1 in p° and <j) production [TT1IT21IT3]. 

We first consider the semi-inclusive production of pions. Depending on the pion charge, exclusive 
channels contributing here are direct production ep — > eir + n and the production and decay of p and 
If* (892). The p decays to almost 100% as p° -> vr+vr- and p + -> vr+vr , and the iP(892) decays to 
almost 100% into Kir, with branching fractions 

B(K* + ^K + 7T°) = i, B(K*+ -> K°7T+) = §, 

B(K* ^K + 7r~) = |, B(K*°^K 7t°) =i (31) 
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Figure 11: Exclusive contributions to the 7r + electroproduction cross section at Q 2 = 2.5 GeV 2 , 
obtained from our leading- twist calculation of a^. This and the following plots are for a lepton beam 
energy of 27.5 GeV in the target rest frame. Left: direct exclusive production, and contributions from 
K* + and p + decay. Right: sum of contributions from K* + and p + decay compared with contributions 
from oj and p° decay. 



following from isospin symmetry. Note that in quark fragmentation one has <r(7r°) = 2 [c( 7r+ )+0"( 7r 
which follows directly from the isospin relations between the pion fragmentation functions. This 
relation also holds for the contributions from K* decay, but it is strongly violated for p decay. For the 
p° this effect was investigated in Ref. [S] in connection with the separation of u and d distributions 
in the proton using semi- inclusive DIS. 

In Fig. we show the result of our leading- twist calculation from Sects. and |S] for the ep cross 
section of tt + production. We give all ep cross section for a lepton beam energy of 27.5 GeV in the 
target rest frame, corresponding to the HERMES experiment, and recall that all our exclusive cross 
sections are calculated with an upper cutoff of 1 GeV 2 on We see that the p° channel is clearly 
dominating. The u, which decays to almost 100% into ir + ir~ir°, is much less prominent. According 
to our discussion in Sects. 0] and H3 one expects that the contribution from ep —>■ en°p to tt° production 
is smaller than in the case of direct tt + production. The same holds for the production and decay 
of r] and 7/, which have several three-body decays contributing to all three pion channels. As we 
argued in Sect.[5]the production of /2(1270), which predominantly decays into ir + -K~ and 7r°7r°, may 
contribute at a similar or lower level as p + decays. In Fig. ^] we show the z spectrum arising from 
different vector meson decays. Whereas p decays contribute in almost the entire z range, pions from 
K* decays are limited to z values below 0.7. Note that due to charge conjugation invariance the z 
spectrum from p° decays is identical for the 7r + and ir~ , and by isospin invariance the same holds for 
the 7r + and 7r° spectra from the decay of p + . To illustrate the dependence of the z spectrum on the 
ratio R of longitudinal and transverse meson production cross sections we have taken values which 
correspond to the range measured in p° and 4> production at Q 2 = 2.5 GeV 2 fTTl IT21 IT31 HI] . 

We now compare the contribution from exclusive channels with the semi-inclusive cross section 
obtained from the leading-order expression (|24| for quark fragmentation. We use the LO parton 
densities from MRST2001 and the LO fragmentation functions of Kretzer j^Hj, both at a scale Q 2 = 
2.5 GeV 2 . Let us first take a look at the high-z tail of the spectrum, where direct exclusive production 
contributes. Integrating the semi-inclusive cross section for 7r + production for z > 0.9 at Q 2 = 
2.5 GeV 2 and xb = 0.1, we obtain da / (dQ 2 dxs) = 0.19nbGeV -2 from (|24[). This number should be 
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Figure 12: Contributions to the semi-inclusive pion electroproduction cross section from decays 
of exclusively produced vector mesons, for Q 2 = 2.5 GeV 2 and x B = 0.1. Shown are the results 
corresponding to our leading-twist calculation of <jl for vector meson production, and a value of 
R = olIot chosen as indicated in the figure. Top left: 7r + from p° decay. Top right: ir + from p + 
decay. Bottom left: 7r + from K* + decay. Bottom right: 7r~ from K*° decay. The contribution to ir° 
production from K* + and K*° decays is given by the average of the corresponding curves in the two 
lower plots, according to (pTTj) . 
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Figure 13: The cross section for semi-inclusive electroproduction of tt + (left) and 7r~ (right), as a 
function of z, at xb = 0.1 and Q 2 = 2.5 GeV 2 . The cross sections from quark fragmentation were 
calculated using the LO fragmentation functions of Kretzer |58j and the MRST 2001 LO parton 
distributions. The contributions from vector meson decays were obtained by adjusting our leading- 
twist results for the vector meson production cross sections, as explained in the text. The value taken 
for the p° cross section matches the HERMES measurement jlUj . 

understood as a naive extrapolation: the fragmentation functions are not well known for z close to 
1, and in the above kinematics z > 0.9 corresponds to an invariant mass mx < 1.84 GeV according 
to (|25[). where leading-twist fragmentation can be just marginally valid. Our leading-twist result 
for ep — ► e7r + n gives da / (dQ 2 dxs) = 0.045 nb GeV -2 , which according to our comparison with 
preliminary HERMES data (Sect. |SJ) undershoots the actual cross section by a factor of about 0.4. 
We thus find that for z > 0.9 direct exclusive pion production may be a substantial part of the 
semi-inclusive cross section, but cannot be more quantitative given the uncertainties just discussed. 

For a realistic estimate of exclusive vector meson production we divide our leading-twist results for 
<jl by a factor 7, except for (f> production. According to Table El this brings us close to the HERMES 
measurement for p° production at Q 2 = 2.3 GeV 2 and xb = 0.1, and according to our arguments in 
Sect. IHlit should give a reasonable estimate for the other channels. In other words, we assume that 
the ratio of vector meson channels is sufficiently well described by our leading-twist calculation. Only 
for <j) production do we divide our leading-twist results by a different factor, namely by 15, following 
our comparison in Table [21 with preliminary HERMES data in this channel. One might argue that for 
the production of a K* , which has one light and one strange quark, power corrections are between 
those for the p° and for the <fi, but we refrain from such refinements here. Possible changes by a factor 
of 2 would in fact not change our conclusions regarding the role of K* decays. For a prediction of ctt 
we divide gl obtained as just described by the value R = ctl/&t = 1-2 from preliminary HERMES 
data for p° production in the relevant kinematics [111 112] . except for the cp channel, where instead we 
take R = 0.8 from the parameterization of preliminary HERMES data in jll| I13j. Variation of R as 
shown in Fig. ^1 would not affect the conclusions we shall draw. 

Integrating the p° decay contribution to 7r + production for z > 0.9 at Q 2 = 2.5 GeV 2 and xb = 
0.1, we find da/{dQ 2 dx B ) = 0.13nbGeV" 2 , which is surprisingly close to the extrapolation of the 
fragmentation result given above. The fragmentation formula for 7r~ production gives da/(dQ 2 dxs) = 
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Figure 14: Same as Fig. 111! but for electroproduction of K + . 



0.07 nb GeV~ 2 when integrated over z > 0.9, so that in this case the p° contribution slightly overshoots 
the naive fragmentation result. In Fig. we show the z spectrum of semi-inclusive tt + and tt~ 
production, comparing the fragmentation result with the contributions from vector meson decays. 
Following our above discussion we do not show the cross section from fragmentation for z above 0.9. 
We see that p° production gives a sizable contribution to semi-inclusive production for z greater than 
about 0.8. According to our estimate, p + production is suppressed relative to p° by two orders of 
magnitude and cannot compete with the cross section from quark fragmentation even at large z. 
The K* decay contribution is somewhat larger in size but limited to z below 0.7. The fragmentation 
result for 7r° production is just the average of 7r + and tt~ because of isospin invariance. With p° decay 
being absent and the contributions from other vector mesons being comparatively small, we find no 
exclusive channel that is prominent in semi-inclusive 7r° production for the kinematics discussed here. 
We expect direct exclusive production ep — > e7r°p to be much less important at high z than in the case 
of 7r + production, to the extent that power corrections enhance the ir + but suppress the 7r° compared 
with the leading approximation in 1/Q 2 - Given the relative size of cross sections in Fig. is clear 
that pions from uj production are significantly smaller than the fragmentation result for all z, and we 
shall not analyze the kinematics of the corresponding three-body decay here. 

Turning to semi-inclusive K + and K~ production, we show in Fig. El the contributions of the 
relevant exclusive channels to the ep cross section, obtained from our leading-twist calculation in 
Sect. El The production of <f>, which decays to approximately 50% into K + K~, is clearly dominant 
for K + production, and it is the only channel contributing to K~ production. As is seen in the 
z-spectra of Fig. 1151 it is however only the K + from K* decays that extends to z values above 0.65. 

Integrating the leading-order fragmentation formula (|24|) for K + production for z > 0.9 at Q 2 = 
2.5 GeV 2 and xb = 0.1, we obtain da / (dQ 2 dxs) = 0.048 nb GeV -2 , which is to be understood as an 
extrapolation as in the pion case discussed above. Our leading-twist estimate for direct exclusive K + 
production in Sect. 21 gives da/(dQ 2 dx B ) = 0.016nb/ GeV 2 . Following our discussion in Sect. one 
expects that power corrections will lead to weaker enhancement than in the case of ep — > e7r + n (or 
possibly even to a suppression), because Ti. is more important in the leading- twist cross section for 
K + production than 8. Nevertheless, the above numbers suggest that direct exclusive K + production 
may be of some significance at the high-z end of the spectrum. 

If we integrate the K + spectrum from K* + and K*° decays for z > 0.9, dividing our leading- 
twist result by 7 and accounting for the transverse cross section as described above, we obtain 
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Figure 15: Same as Fig. 1121 but for the contribution from K* + and K*° decays to K + production 
(left) and for the contribution of <fi decay to K + or K~ production (right). 



da I (dQ 2 cIxb) = 0.0021 nb GeV~ 2 , which is well below our extrapolation from leading-twist frag- 
mentation. In Fig. E3 we compare the fragmentation result for semi-inclusive K + and K~ production 
with the individual contributions from K* and 4> decays. We conclude that, even within the uncer- 
tainties of our estimates, contributions from K* production are only a fraction of the fragmentation 
result at any z, and that <p production, despite its larger cross section, is always well below the 
semi-inclusive cross section. Our finding concerning the (ft contribution agrees with a recent study 
of measured kaon multiplicities in [H]. On one hand, kaon production by quark fragmentation is less 
suppressed compared with pion production than is exclusive <j) production compared with production 
of p°, and on the other hand (j) decays only contribute in a z-range where the fragmentation functions 
are still large. Apart possibly from direct K + production at z close to 1, we thus find no exclusive 
channel dominating K + or K~ production in typical HERMES kinematics. 

So far we have compared exclusive channels with quark fragmentation at xb = 0.1 and Q 2 = 
2.5 GeV 2 . We have also performed the comparison of Figs. El and El for xb = 0.3 at the same Q 2 . 
For the vector meson cross sections we used the same values of R and the same correction factors 
as for xb = 0.1, dividing our leading-twist cross sections by 7 for all vector mesons except the </>, 
where we divide by 15. In doing this, we assume that the leading-twist approximation gives a realistic 
description of the xb dependence in this region. We find that our qualitative conclusions do not 
change when going to the larger value of xb- 

A comment is in order concerning the treatment of exclusive channels in the analysis of semi- 
inclusive DIS data when extracting quark fragmentation or distribution functions. It is by no means 
clear that by subtracting contributions of individual exclusive channels from the total yield one 
obtains an observable more suitable for comparison with the quark fragmentation formulae. In fact, 
the derivation of factorization theorems relies on the sum over all channels X in (|22j) to be complete. 
At sufficiently large Q 2 , each individual exclusive channel by itself is a power correction which may 
or may not be included in the leading-twist analysis. This situation is similar to the one with the 
contribution from individual nucleon resonances to the cross section of inclusive DIS. A way to think 
about the relation of exclusive channels to the leading-twist cross section is quark-hadron duality. 
It remains a challenge to formulate the problem of quark-hadron duality for semi-inclusive DIS in a 
quantitative fashion. 

Symmetry properties like a(ir°) = ^[a(7r + ) + a(TT~)] can emerge from summing over many channels 
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Figure 16: Same as Fig. 1131 but for semi- inclusive production of K + (left) and K (right). As 
discussed in the text, the <fi production cross section used in this plot matches the preliminary results 
of a HERMES measurement PUE3]. 

which do not fulfill this relation individually. If however a single channel like p° production domi- 
nates the semi-inclusive cross section, it clearly becomes more and more difficult for the remaining 
channels to compensate the missing symmetry between charged and neutral pion production. In such 
a situation, parton-hadron duality must cease to work. The outcome of our study is that indeed the 
only channel whose contribution can become dangerously large in HERMES kinematics, is the p° 
contribution to ir + and tt~ production. 

8 Summary 

We have evaluated the cross section for a variety of exclusive meson production channels for moderate 
to large xb at leading order in 1/Q 2 and in a s . Cross sections change significantly when varying the 
nonperturbative input, generalized parton distributions and meson distribution amplitudes, within 
plausible limits of current model building. On one hand this implies an uncertainty in predicting 
these cross sections, but on the other hand it implies that their measurement can ultimately help to 
constrain the nonperturbative functions, provided theoretical control over corrections to the leading- 
order formulae. We find the largest cross section uncertainties for p°, u> and (p production, which 
is sensitive to the generalized gluon distribution over a large range of xb, reflecting the current 
uncertainty of the unpolarized gluon density at low scales. A strong dependence on the factorization 
scale in these channels underlines the need for analysis at next-to- leading order in a s . Comparing 
our leading-twist cross section with experimental data, we confirm that for Q 2 of a few GeV 2 power 
corrections are substantial. In particular, the suppression of vector meson cross sections we find 
is consistent with what has been estimated in the literature from the effects of parton transverse 
momentum in the hard-scattering subprocess. A consistent description of such effects together with 
next-to-leading order corrections in a s remains a challenge for theory. As is seen for the ratio of <p 
and p° production, the most serious theoretical uncertainties cancel however in cross section ratios for 
sufficiently similar channels (the main distinction being between channels with and without i-channel 
pion exchange). 
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Rescaling our leading-twist cross sections such as to be consistent with experimental data for p 
and (ft production, we have compared their contribution to semi-inclusive pion or kaon production 
with the result of leading-twist quark fragmentation, focusing on the typical kinematics of HERMES 
measurements, where Q 2 ~ 2.5 GeV 2 and xb ~ 0.1. Within large uncertainties, direct exclusive pro- 
duction of 7r + and possibly K + appears to be comparable with the fragmentation result extrapolated 
to the bin 0.9 < z < 1. Through their decays, exclusively produced p, (ft and K* contribute in a wide 
range of z. Pions from K* decay and kaons from (ft decay are however limited to z below 0.7. With 
this and the relative size of cross sections, our estimates indicate that in typical HERMES kinematics 
the only exclusive channel whose cross section can compete with quark fragmentation is the p°. The 
p° saturates the quark fragmentation result for semi-inclusive tt + and ir~ production at large z. Since 
the p° does not contribute to ir° and to kaon production, there is no corresponding "dangerous" vector 
channel in these cases. 
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A Integrals over GPDs within the double distribution model 

The t independent functions in the ansatz ()14|) for quark and gluon GPDs are modeled as 



H«(x,0 



dp 



d(3 



1+1/31 



da 5{x-f3- £a) h{(3, a) 9{(5) q(j3) - 0(-p) q(-0) 



1+1/31 



da 5{x- (3- fa) h(/3, a) (3 0(J3) g((3) - 6{-(3) g{-[3) 



-i 



dfl 



-i+ 



da 5{x-f3- £a) h(p, a) 9{(3) Aq(j3) + 9{-p) Aq(-0) 



(32) 



where 6 denotes the usual step function, q, q, Aq, Aq the unpolarized and polarized quark and 
antiquark distributions, and g the unpolarized gluon distribution. The profile function 
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(33) 



depends on a parameter b, which we chose to be either b = 1 or b = 2 in this work. 
For meson production amplitudes we need the integrals 
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where we used the definitions H q (x,^) = —H q (—x,^) and H q (x,t;) = H q (—x,t;) together with the 
fact that these functions are even in f . The required integral for gluons can be brought into the same 
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Figure 17: The functions appearing in the integrals ()36|) and (j37|) for the meson production amplitude 
for £ = 0.053 (top) and £ = 0.18 (bottom), which respectively correspond to xb = 0.1 and xb = 0.3. 
Note the different ranges for (5 in the plots. 



form as I q {Cj by rewriting 



ni)= i\ x H '^\ 1 . -U 1 , (35) 



i x £ — x — £ l £ — f — ie 

where we used that H 9 (x, £) is even in x. The imaginary parts of these integrals are readily converted 
into integrals over (3, with 

r— 

Iml^O = dpi(j3,t)q(l3), 
US,) = ^(26 + 2) (l-a b 1 (JL_3) b B b (36) 

'W'V 2 2& + i r 2 (6+i) e h+1 (i-/?) 26+i vi+e ') 1 { ' 

for £ > 0. The function /(/?, £) vanishes at the endpoints of the integration region, which in particular 
ensures the convergence of the integral at /3 = for common parameterizations of quark densities. 
To obtain the analogous expressions for fl, I q , I 9 and I 9 one has to replace q{(5) with q(j3), Aq(/3), 
Aq(f3) and (3g(f3), respectively. 

The real parts of the amplitudes involve principal value integrals, whose numerical evaluation 
requires some care, especially for small £. For our choices of profile parameters 6=1 and 6 = 2 one 
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can explicitly perform the a integral after inserting (|32|) into (|34j) and (jHSJ- The result is 
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For both 6 = 1 and 6 = 2, the function i?(/5, is continuous in the full interval of integration, with 
finite limits at /5 = and (3 = 1. If £ > it is positive for (3 < £ and negative for (3 > £, and if 
£ < it is negative in the entire interval. Convergence of the integral for polarized quark distributions 
requires that Aq(f3) and Aq(f3) have integrable singularities at (3 = 0, which is the case for the parton 
densities we use in this study. The unpolarized quark distributions have a steeper behavior at small 
(3, but since i?(/5,£) + i?(/5, — £) ~ /5 for /5 — ► it is sufficient to have integrable singularities for 
?(/?) - ?(/5) and for [3q{(3). 

In Fig. El we illustrate the behavior of the functions multiplying the parton distributions in the 
integrals (|36j) and l|37|l. The imaginary part of the amplitude involves momentum fractions in the 
parton densities between and 2£/(l + £) = xb, with a maximum of the shape function I{(3, £) for 
/3 around £. In contrast, the real part is sensitive to higher momentum fractions, with a partial 
cancellation from values above and below £. One also clearly sees the stronger sensitivity to small (3 if 
6=1. Note that the functions shown in the figure will be multiplied in the amplitude with functions 
showing a strong rise towards (3 = 0. 



B Distribution of pions or kaons from vector meson decay 

In this appendix we discuss the decay of a vector meson into two pseudoscalar mesons and derive the 
z distribution given in (|27|l . Consider the contribution of ep — ► VB with subsequent decay V — > P1P2 
to semi-inclusive production ep — > P\ + X. A useful set of variables to describe the decay of the 
vector meson are the polar and azimuthal angles 9 and 93 of P\ in the vector meson center-of-mass, 
as shown in Fig. 1181 The distribution in these angles is connected in a straightforward way with the 
spin density matrix of the produced vector meson |57| . and the phase space element has a factorized 
form in the variables Q 2 , xb, t and 9, <p. The variable z used for semi-inclusive production of P\ is 
then given by 

z = a + 6 cos 9 + c sin 9 cos ip (39) 



29 




V rest frame 



Figure 18: Kinematic variables for 7*(g) + p(p) — > V(q') + B(p') followed by the decay V — > P1P2, 
shown in the 7* p center-of-mass and in the rest frame of V. Here 9 and </? are the spherical coordinates 
of the momentum of Pi in the depicted coordinate system. 
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where we abbreviated 
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The energy and momentum Ep\ and \Qpi\ of Pi in the rest frame of V have already been given in 
(j30j) . and At is the transverse momentum of the scattered baryon with respect to the initial proton 
in the 7*p center-of-mass (see Fig. I18|) . The approximate expressions in (|4U|) are valid up to relative 
corrections of order x B m 2 /Q 2 and x B A^/Q 2 , and to the same accuracy one has A 2 ^ = (1 — x B )(tQ — t). 
Changing variables from 9 to z gives for the cross section 
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(42) 



In Bjorken kinematics one has c <C 6 and can replace the Jacobian in (|42j) by 1/6 (except in the small 
region where sin 9 ~ c/b, which is not relevant for our purposes). Neglecting Ap we get z = a + b cos 9 
with a and b evaluated at Ap = 0, and integration over t and (p gives 
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in terms of the cross sections for the production of longitudinally or transversely polarized vector 
mesons. Using s-channel helicity conservation, which is experimentally seen to hold at the few 10% 
level in p° and <fi production |11 | I12 | IT3*]. these cross sections respectively correspond to the production 
from longitudinally or transversely polarized virtual photons, and we finally obtain ()27|) . In our 
numerical applications we have used the exact expressions from lf4"0"|) and (|4*Tj) at At = 0, and thus 
in particular neglected c. Since the integrated cross sections are dominated by small At, this should 
be a very good approximation for the values of Q 2 and xb we focus on in the present study. The 
inclusion of finite effects in the kinematics would considerably complicate any analysis. 
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